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Abstract 

Land planarians are an interesting group of free-living flatworms that can be useful as bloindicators because of their high 
sensitivity to environmental changes and low dispersal capacity. In this study, we describe and compare assemblages of 
land planarians from areas with different conservation degrees of the Interior Atlantic Forest (Misiones, Argentina), and 
assess factors that could be related to their abundance and richness. Eight sites were tracl<ed in search of land planarians in 
Reserva de Vida Silvestre Urugua-f (RVSU) and Campo Anexo Manuel Belgrano (CAMB). Diurnal and nocturnal surveys were 
performed in each site along nine sampling campaigns. We collected 237 individuals belonging to 18 species of the 
subfamily Geoplaninae. All sites were dominated by Geoplana sp. 1 and Pasipha hauseri. The richness estimators showed 
that there would be more species in RVSU than in CAMB. The abundance and richness of land planarians was high during 
the night and after rainfalls, suggesting an increased activity of flatworms under such conditions. The abundance and 
richness of land planarians were also related to the conservation condition of the sites. Disturbed sites showed less 
abundance and richness, and were segregated from non-disturbed ones by nmMDS analysis. Beta diversity between sites 
was higher than expected, indicating that the species turnover between sites contributed more to the total richness 
(gamma diversity) than the alpha diversity. 
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Introduction 

Land planarians (Platyhelminthes: Geoplanidae) successfully 
colonized the terrestrial environment mUlions of years ago [1]. 
However, they have not been able to develop mechanisms for 
water conservation and are thus unable to withstand desiccation 
[2]. In addition, land planarians have low vagUity and cannot 
endure long periods of immersion in water, so they can be 
considered stenoic organisms regarding their habitat requirements, 
i.e., they are very sensitive to the moisture conditions of the 
environment [2] . Land planarians are 'top predators' of the soil 
fauna [3,4]. They can feed on a wide range of soil invertebrates, 
mainly earthworms, snails, lugs, leeches, insects, isopods, and 
arachnids [5—10]. Therefore, they may be good environmental 
indicators, particularly in tropical and subtropical rainforests, 
where they are abundant [3]. 

Land planarians exhibit the highest diversity in the Atlantic 
Forest, with about 180 species described for its Brazihan portion 
[11]. The Atlantic Forest is one of the world's 25 recognized 
biodiversity hotspots [12]. It is a complex of ecosystems which 
extends along the Atlantic coast of Brazil and inland as far as 
eastern Paraguay and north-eastern Argentina, constituting the 
Interior Adantic Forest. The original coverage of the Atlantic 
Forest remains in small fragments and under some kind of 
conservation status [13,14]. In Argentina, there are stiU large 



extensions of the original Adantic Forest, mainly under some kind 
of legal protection. This ecosystem has been suffering from human 
impact that has modified the original landscape, due to 
uncontrolled deforestation (which has been the primary cause of 
forest degradation), the burning of the land to prepare it for 
farming or grazing, and the introduction of exotic species with 
commercial purposes [13]. 

Several researches have study the diversity of land planarians in 
the Atlantic Forest of Brazil [15-18], and some have compared 
different assemblages in man-disturbed areas [5,19]. At present, 
the diversity of land planarians in the Argentine portion of the 
Adantic Forest is unknown. Recendy, we have started to describe 
the diversity of these flatworms in this region [20,21]. In order to 
improve our knowledge about land planarians, in this work we 
describe and compare assemblages from areas of the Interior 
Adantic Forest of Argentina with different conservation degrees, 
and assess factors that could be related to their abundance and 
richness. Additionally, we determine the contribution of species 
richness (alpha) and species turnover (beta) to the gamma diversity. 

Methods 

Ethics Statement 

The research has been conducted according to the Argentine 
law. 
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Study area and sampling design 

Our study took place in the southern portion of the Interior 
Atlantic Forest. In Argentina, this region c:overs about 25,700 km^ 
and is part of the Paranaense Subregion (Neotropical Region), 
where two biogeographical provinces are recognized: the Para- 
naense Forest and the Moist Forest with the coniferous tree 
Araucaria angustifolia [22]. These provinces are characterized as 
semi-deciduous subtropical rainforest. The climate is warm and 
humid, with an annual temperature of 16-22°C and a total annual 
precipitation of 1,600-2,000 mm [23]. 

The surveys were performed in two reserves of Misiones 
province, separated by about 40 km: Reserva de Vida SUvestre 
Urugua-i (RVSU) (25° 59' S, 54° 05' W) and Campo Anexo 
Manuel Belgrano (CAMB) (26° 02' S, 53° 47' W), each 
representing the Paranaense Forest and Moist Forest with A. 
angustifolia, respectively (Figure 1). The reserves are differentiated 
by conserved surface, altitude, management degree, and vegeta- 
tion type. RVSU is a private natural reserve that covers 3,423 ha 
at ~200 m a.s.l. It was created in 1997 and previously used for 
selective logging until the 1970s. This reserve, now under strict 
protection, is part of one of the largest corridors of continuous 
original rainforest in the southern portion of the Atlantic Forest, a 
'green block' of almost 6,000 km^ [24]. It is characterized by 
diversified forests, although trees of Balfourodendron riedeliamm and 
Nectmdra spp. dominate plant formations. CAMB is a governmen- 
tal forest reserve that covers 2,136 ha at ~600 m a.s.l. It was 
created in 1948 to protect native and planted populations of ^. 
angustifolia. This rainforest is also characterized by an undergrowth 
of tree ferns (Alsophyla sp., Dicksonia sp., Trichipteris sp.) [24]. 
However, in CAMB there are also plantations with exotic conifers 
{Pinus taeda). Therefore, this reserve is a mosaic of preserved and 
disturbed areas, isolated from other protected areas and 
surrounded mainly by small farms (Figure 1). 

In each reserve, four sites were tracked in search of land 
planarians. All sites were selected for their accessibihty, using 
existing paths. In RVSU, the four sites (U1-U4) were located in 
undisturbed areas with a condition of native vegetation and with 
the same physiognomy. In CAMB, the four sites (B1-B4) were 
situated in heterogeneous land use areas. Sites Bl and B3 were 
situated in undisturbed areas with A. angustifolia, while B2 and B4 
were located in disturbed areas with exotic vegetation [P. taeda) 
(Figure 1). 

In both reserves, land planarians were collected with the same 
sampling procedure along nine campaigns, between 2008 and 
2010. In each site, samplings were performed by one person 
walking along the paths (1,000-1,500 m long) for 2 h during the 
day and another 2 h during the night. Therefore, 36 hours of 
sampling (2 h diurnal x2 h nocturnal x9 campaigns) were carried 
out in each site, totaling 144 hours of sampling effort for each 
reserve. During the day, the collector searched for planarians 
beneath and inside fallen logs, leaf litter and stones, whereas 
during the night, when planarians are more active, the collector 
performed the direct observation of the soU by means of a 
headlamp. Land planarians were manually collected and placed in 
plastic recipients with humid leaf litter, to avoid dehydration stress. 

Environmental data were obtained from climatological stations 
near the rescr\'i's. The mean temperature of each sampling day 
and the cumulative rainfall (sum of millimeters of the ten days 
previous to each sampling day) were recorded. 

Morphological analysis 

The external features of each specimen were observed in live. 
Planarians were then killed by throwing boiling water on them to 
avoid distortions of their tissues and then fixed in 10% 



formaldehyde and preserved in 70'/o ethanol. Fragments of 
different body regions (anterior region, pre-pharyngeal region, 
pharynx and copulatory apparatus) were dehydrated and embed- 
ded in Paraplast, serially sectioned at 6-10 [im thick with a 
microtome, and stained with hematoxylin-eosin and the Masson's 
trichrome method [25]. The histological preparations were 
observed by optical microscope and the copulatory apparatus 
reconstructed for identification purposes. The specimens were 
studied by two specialists on free living flatworms (LN and FB). 
Some specimens were identified at species level and others as 
morphospecies because they are in description process. 

The material studied was deposited in the Invertebrate 
Collection at Museo de La Plata (MLP), Argentina. 

Data analysis 

Species richness is the simplest way to describe the diversity of a 
community and to make comparisons [26]. We disaggregated 
species richness in (a) punctual alpha diversity, for the number of 
land planarian species recorded in each site, and (b) cumulative 
alpha diversity, for the total number of species in each reserve. We 
constructed species accumulation curves using sample-based 
rarefaction, with the expected richness function Mau Tau (with 
95'!^) confidence intervals), to compare cumulative alpha diversity 
curves between reserves [27]. We used the species accumulation 
curves with the number of individuals instead of the number of 
samples to avoid biases in comparison due to differences in the 
abundance of land planarians [28]. Since all inventories have 
unrecorded species [29], to analyze and compare the completeness 
of the species inventory in each reserve, we tested the performance 
of eight richness estimators, based on abundance (Chao 1, ACE) 
and incidence (Chao 2, ICE, jacknife 1, jacknife 2, bootstrap, and 
Michaelis-Menten), and the number of singletons and doubletons. 
Species accumulation curves and richness estimators were 
computed using Estimates v.9.0 [30], performing 100 randomiza- 
tions in each analysis. To analyze the dissimilarity among reserves, 
we performed the calculation as 1 - Bray-Curtis index, using 
Chao's procedure, with Estimates v.9.0, because it compensates 
the effect of unseen shared species [31]. 

Abundance and species richness are interesting parameters to 
assess the structure of assemblages, since they are simple and fast 
to measure. Furthermore, these parameters can change according 
to environmental conditions. Besides the effects of mean temper- 
ature and cumulative rainfall on abundance and richness of land 
planarians, we also evaluated the relationship of these parameters 
with the conservation condition of each site. Therefore, we 
performed a multiple regression analysis (GLM), including data of 
both reserves (2 reserves x4 sites x9 campaigns), to test the effects 
of three independent variables: mean temperature, cumulative 
rainfall (continuous predictor variables), and conservation condi- 
tion (categorical predictor variable: disturbed/undisturbed) on two 
dependent variables (abundance and richness). 

To assess whether we had found more land planarians in 
diurnal or nocturnal surveys, we used a non-parametric WUcoxon 
test, for dependent samples. We compared the abundance of 
planarians collected in 36 diurnal and nocturnal samplings (4 sites 
x9 campaigns) in each reserve. After this comparison, we pooled 
diurnal and nocturnal sampling data, summarizing nine replicates 
(campaigns) per site. This procedure ensured that the actual 
richness and abundance of each site was well estimated. 

We then used a two-way ANOVA to tested the effect of the 
reserves (fixed, orthogonal factor, with two levels: RVSU and 
CAMB) and sites (random, nested factor, with U1-U4 nested in 
RVSU level, and B1-B4 nested in CAMB level) on the abundance 
and richness of land planarians. Data were log-transformed to 
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Figure 1. Study area, in tKie Argentine portion of tlie Interior Atlantic Forest (Misiones province). 1, Reserva de Vida Silvestre Urugua-f 
(RVSU), with four undisturbed sampling sites (U1-U4), and 2, Campo Anexo Manuel Belgrano (CAMB), with two undisturbed sampling sites (B1 and 
B3), and two disturbed sampling sites (B2 and B4). Maps modified from Di Bitetti et al. (2003), Galindo-Leal and Camara (2003), and Google Earth. 
doi:l 0.1 371/journal.pone.009051 3.g001 



fulfill ANOVA assumptions. Post-hoc Student-Newman-Keuls 
(SNK) tests were applied for multiple comparisons. 

We buUt rank-abundance (dominance-diversity) curves as 
indicators of the structure of the planarian assemblages of each 
site in both reserves. The relative abundance of each taxon on a 
logarithmic scale (logio) was plotted against the rank order of the 
taxa, from the most to the least abundant. These curves constitute 
a useful tool to visualize some aspects of the assemblages such as 
species richness (number of points), evenness (slope), number of 
rare species (tail of curves) and relative abundance of each species 
(order of the species in curves) [32]. The composition of land 
planarian assemblages was compared among sites of RVSLI and 
CAMB by the non-metric multidimensional scaling (iimMDS), 
using the Bray-Curtis coefficient, to evaluate the similarities based 
on the abundance matrix (log-transformed). 

Additionally, we determined the degree of contribution of 
species richness (alpha) and species turnover (beta) between sites to 
landscape scale diversity (gamma diversity). We adopted the 



framework of Jost [33,34] to calctilate alpha (diversity within 
samples) and beta (diversity among samples), in which gamma 
diversity (total amount of diversity) is partitioned according to the 
formula H^, *Hp = H,^ (H: Shannon- Wiener entropy). Entropies, 
like the Shannon- Wiener index, are not themselves diversities, and 
their use may obscure differences in diversity because indices differ 
only by small magnitudes [35] . Therefore, we used a transforma- 
tion that allows an intuitive interpretation of species diversity by 
using the effective number of species, named by Jost [33], as a 
measure of "true diversity". We computed the diversity with q = 0, 
where diversity is completely insensitive to the abundances of 
species and the value obtained is thus equivalent to species richness 
("Z) = Sobs), and g= I, the exponential Shannon entropy, where 
each species is included with a weight proportional to its 
abundance ('D=exp H') [36]. We used Partition v.3.0 [37], 
using 1000 randomizations (Monte Carlo resampling method) to 
derive the expected values of alpha and beta diversity that would 
be obtained if individuals or samples were randomly distributed. 
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Table 1. Abundance and 
the species in Figure 4. 


richness (alpha diversity) of land planarian 


species collected 


in each sampling site. Letters (A-R) identify 




Reserves 


RVSU 




CAMB 










Taxa/Sampling sites 


U1 


U2 U3 U4 


Total B1 


B2 B3 


B4 


Total 




Choeradoplana crassiphalta 


A 3 


5 


8 


2 


1 


3 


11 


Geoplana sp. 1 


B 18 


38 7 12 


75 10 


1 13 


3 


27 


102 


Geoplana sp. 2 


C 1 


6 


7 4 


7 


1 


12 


19 


Geoplana sp. 3 


D 






2 




2 


2 


Geoplana sp. 4 


E 


1 


1 


2 




2 


3 


Geoplana sp. 5 


F 


2 


2 








2 


Geoplana sp. 6 


G 




1 


1 




2 


2 


Geoplana sp. 7 


H 






1 




1 


1 


Gigantea sp. 1 1 - 1 - - 1 _____ ■] 


Matuxia cf. matuta 


J 




1 






1 


1 


Pasipba hauseri 


K 31 


5 2 4 


42 11 


9 10 




30 


72 


Pasipha sp. 1 


L 1 




1 








1 


Pasipha sp. 2 


M 


3 


3 








3 


Pasipba sp. 3 


N 2 




2 








2 


Pasipha sp. 4 


0 






1 




1 


1 


Pasipba sp. 5 


P 






1 




1 


1 


Supramontana argentina 


Q 1 


1 - 1 


3 


2 




2 


5 


Xerapoa cf. pseudorhynchodemus 


R 2 


1 2 


5 2 




1 


3 


8 


Total abundance 


59 


46 12 33 


150 29 


10 42 


6 


87 


237 


Punctual alpha diversity 


8 


5 4 7 


6 


2 11 


4 






Cumulative alpha diversity 


12 




13 











doi:l 0.1 371/journal.pone.009051 3.t001 

Results 

We found 237 individuals (150 in RVSU and 87 in CAMB), 
representing 1 8 land planarian species distributed in six genera of 
the subfamily Geoplaninae. We found 12 species in RVSU and 13 
species in CAMB (Table 1). For cumulative alpha diversity, the 



species rarefaction curves were not different between reserves 
(Figwe 2). However, singleton and doubleton curves showed 
different trends. In CAMB but not in RVSU, singleton and 
doubleton curves reached the intersection with each other 
(Figure 2). Among the richness estimators analyzed, bootstrap 




160 



Individuals 



Figure 2. Land planarian species rarefaction curves, singleton and doubleton curves for reserves. Solid lines show the rarefaction curves 
(bars delineate 95% confidence intervals). Dashed and dotted lines represent the singleton and doubleton curves, respectively. RVSU (squares) and 
CAMB (circles). 

doi:10.1371/journal.pone.0090513.g002 
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Table 2. Land planarian richness, singletons and doubletons 
observed in each reserve. 





RVSU 


CAMB 


Cumulative alpha diversity 


12 


1 3 




Number of singletons 


4 






Number of doubletons 


2 


4 




Chao 1 


17 71% 


15 


87% 


ACE 


1 6 75% 


16 


81% 


Chao 2 


28 43% 


16 


81% 


ICE 


1 9 63% 


20 


65% 


Jacknife 1 


17 71% 


17 


76% 


Jacknife 2 


21 57% 


19 


68% 


Bootstrap 


1 5 80% 


15 


87% 


Michaelis IVlenten 


1 9 63% 


19 


68% 


Average of inventory completeness 


65% 




77% 



Number of species expected and percentages of inventory completeness 
according to different richness estimators. 
doi:10.1371/journal.pone.0090513.t002 



was the most consei'vative one, while Chao 2 estimated the highest 
number of species in RVSU (Table 2). In CAMB, Chao 1 and 
bootstrap estimated the lowest number of species, similar to that 
recorded, whereas ICE estimated the highest number of species. 
On average, taking into account the eight estimators, we reached 
65% of completeness of the inventory of species for RVSU and 
77% of that for CAMB (Table 2). Five species were unique to 
RVSU, whereas six species were unique to CAMB, being seven 
species shared by both reserves (Table 1). According to Chao 
estimation, we found a dissimilarity of 38% between reserves 
(Bray-Curtis similarity index = 0.62). 

During the study, the mean temperature varied between 9.6°C 
and 27.9°C, and the cumulative rainfall varied between 6.5 mm 
and 1 30 mm. The cumulative rainfall showed a positive effect on 
the abundance (p = 0.037) and richness (p = 0.020) of land 
planarians, while the temperature was not related. The conserva- 
tion condition also affected the abundance and richness (Table 3). 
In both reserves, we found a greater abundance of land planarians 
during the nocturnal samplings (RVSU: T= 33.5, p<0.0001 and 
CAMB: T= 14.5, p<0.001) (Figure 3). Approximately 90% of the 
individuals collected in RVSU (N=137) and 80% of those 



Table 3. Summary results of a GLM - multiple regression 
analysis, testing the relationship between independent 
variables (continuous: mean temperature, cumulative rainfall; 
and categorical: conservation condition) and two dependent 
variables (abundance and richness of land planarians). 







Dependent variables 






Abundance 




Richness 




Independent variables 


df 


MS 


F 


P 


MS 


F 


P 


Mean temperature 


1 


2.67 


0.184 


0.669 


2.53 


1.123 


0.293 


Cumulative rainfall 


1 


65.15 


4.493 


0.037 


12.75 


5.653 


0.020 


Conservation condition 


1 


110.2 


7.602 


0.007 


22.16 


9.817 


0.002 


Error 


68 


14.5 






2.25 







collected in CAMB (N= 70) were recorded during the nocturnal 
surveys. 

Geoplana sp. 1 and Pasipha hau.seri (Froehlich, 1959) together 
represented 74% of the total collected planarians (Table 1). 
Geoplana sp. 2, Choeradoplana crassiphalla (Negrete & Brusa, 2012) 
and Xerapoa cf pseudorhynchodemus showed an intermediate level of 
abundance, totaling 16% of the land planarians collected. The 
remaining species showed low abundance (Table 1, Figure 4). The 
general abundance of land planarians varied from 12 to 59 
individuals in RVSU, and from 6 to 42 in CAMB. In RVSU, the 
richness ranged from four species in U3 to eight species in Ul. 
Geoplana sp. 1 was dominant in sites U2-U4 and Pasipha hauseri was 
dominant in Ul. Only these two species were common to Ul— U4. 
In CAMB, the number of species varied from two in B2 to 1 1 in 
B3. Only Geoplana sp. 1 was recorded m the four sites, being 
dominant in B3 and B4, and Bl together with P. hauseri, while the 
latter was dominant in B2 (Table 1, Figure 4). 

The reserves were not an important factor on the abundance 
(F= 1.789, p = 0.229) or on the richness (F = 0.789, p = 0.409) of 
land planarians. However, the abundance and the richness were 
sensitive to the different sites (Table 4). The disturbed sites B2 and 
B4 showed a lower abundance and richness than the undisturbed 
sites (U1-U4, Bl and B3) (Figure 5). Based on the composition of 
land planarian assemblages, the nmMDS analysis clustered sites 
with undisturbed condition at 54% similarity (Bray-Curtis 
coefficient). This analysis also segregated B2 and B4, which are 
disturbed sites (Figure 6). 

Partitioning of the gamma diversity revealed that alpha species 
richness and alpha exponential Shannon diversity were not 
dilferent from expected. However, beta diversity between sites 
was higher than expected for both measurements (Table 5). 

Discussion 

We observed no differences in cumulative species richness 
curves between reserves, and rising curves imply incomplete 
inventories. However, the singletons curve in CAMB reached the 
asymptote, tended to zero and intersected the doubletons curve, 
indicating that the species inventory in this reserve would be close 



20 
Q- 18 

E 

« 16H 
ui 

ro 14 



12- 



10 



o 

0) 

o 
c 
ro 

"D 

c 

< 



RVSU 



N=137 



Wilkoxon 

T= 33.5 p<0.001 



N=13 



CAMB 

Wilkoxon 

T= 14.5 p<0.001 



N=70 



N=17 



Diurnal Nocturnal Diurnal 

Sampling period 



Nocturnal 



doi:1 0.1 371/journal.pone.009051 3.t003 



Figure 3. Comparison of land planarian abundance between 36 
pairs of diurnal and nocturnal surveys, in each reserve (RVSU 
and CAMB). Non-parametric Wilcoxon test, for dependent samples 
was used. Mean (dot), median (line inside the box), 25-75% (box), min- 
max (whiskers), and N= total abundance of land planarians. 
doi:10.1371/journal.pone.0090513.g003 
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Figure 4. Rank-abundance curves of land planarian assemblages in each sampling site. Species codes are given In Table 1. Geoplana sp. 1 

(B) and Pasipha hauseri (K) were the most abundant species In all sites. 

dol:10.1371/journal.pone.0090513.g004 



to completeness [38]. The species richness estimators are 
consistent with this assumption. In contrast, in RVSU, the 
singleton and doubleton curves suggest that the species inventory 
is still far from being complete. The species richness estimators for 
RVSU were more variable than those for CAMB. With Chao 1, 
which is sensitive to rare species [39], we reached about 70% of 
the species inventory in RVSU and almost 90% in CAMB. 
Nevertheless, according to Chao 2, which takes into account 
presence/absence data, we reached about 40% of the species 
inventory in RVSU and over 80% in CAMB. Therefore, it would 
be expected that not observed-species are recorded in new surveys. 
Other contributions on land planarian diversity did not reach the 
completeness of the species inventory either. Carbayo et al. [5], 
using Chao 1, and Leal-Zanchet et al. [15], using Chao 2, 
obtained 60% of the species inventory for chflFerent land planarian 
assemblages. Considering that RVSU is part of the largest corridor 
of continuous preserved rainforest, it would be expected that this 
area would have higher species richness than CAMB. However, 
we found similar richness in both reserves, but according to Chao 
2, which is unbiased for small samples [38], 16 species were not 
recorded in RVSU and only three species were not recorded in 
CAMB. 



Studying different land planarian assemblages, Antunes et al. 
[16] recorded a similarity of ~10% in species composition among 
two peri-urban forest areas, distant each other by about 60 
kilometers. Baptista et al. [40] recorded 30% of similarity between 
two land planarian assemblages separated by over 260 kilometers 
in southern Brazil. Compared to these results, the similarity 
estimated between RVSU and CAMB (—60%) is rather high. This 
similarity can be a consequence of a short distance between 
reserves, and due to the fact that RVSU is located on the eastern 
limit of the Paranaense Forest, very close to the Moist Forest of 
Araucaria, which would allow the sharing of species. 

The cumulative rainfall, but not the temperature, affected the 
abundance and richness of flatworms. It is probable that the high 
moisture available everywhere makes unnecessary for land 
planarians to retreat to refuges (e.g. under fallen logs, leaf litter, 
and stones) [41]. Most individuals in both reserves were collected 
during nocturnal samplings. Although the behavior of land 
planarians is stiU poorly understood, our findings may incUcate 
an increased activity of planarians during the night and after 
rainfalls, which would explain the high abundance and richness 
recorded under such conditions. In other studies, in which the 
samplings were carried out only during the day, no relationship 



Table 4. Summary results of a two-way ANOVA model, testing effects of the reserves (fixed; RVSU and CAMB) and sites (random, 
Ul to U4 nested in RVSU and B1 to B4 nested in CAMB) on the abundance and richness of land planarians. 











Abundance 






Richness 












C = 0.218 (NS) 






C = 0.208 (NS) 










df 


MS 


F 


P 


QM 


F 


P 


Reserves 


Fixed 


1 


0.681 


1.789 


0.229 


0.154 


0.789 


0.409 
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Figure 5. Comparison of abundance and richness of land 
planarians between sampling sites (U1-U4 and B1-B4). Means 
and standard errors. Different labels represent significant differences 
(Student-Newman-Keuls; a = 0.05). 
doi:1 0.1 371/journal.pone.009051 3.g005 



was observed between abundance and richness of land planarians 
with environmental conditions [15,42]. 

The pattern of species abundances observed in both reserves, in 
which a few species are ver^' abundant [Geoplana sp. 1 and Pasipha 
hausertj, others are moderately represented, and many are rare 
species, seems to be common for land planarian assemblages, in 
agreement with that previously reported [5,17,18,43]. The high 
abundance of Geoplana sp. 1 and P. hauseri in all sites surveyed can 
be due to the fact that these species are generalists, with the ability 
to colonize undisturbed habitats and habitats with different 
disturbance degrees, being more tolerant or with better capacity 
to thrive in environments under different conditions. Moreover, 
we found these species in pine plantations and gardens in small 
towns close to the reserves. Similarly, other geoplanid species, such 
as Choeradoplana iheringi Graff, 1899, Obama tadislavii (Graff, 1899) 
and Pamba franciscana (Leal-Zanchet & Carbayo, 2001), which are 
the most abundant flatworms recorded in different assemblages of 
southern Brazil [5,18,19,44], have been also found in preserved 
and man-disturbed areas, even close to human setdements, such as 
small gardens in cities and dump deposits [5,44]. 

The abundance and richness of flatworms were related to the 
conservation condition of the sites. The lower abundance and 
richness in the disturbed sites (B2 and B4), is probably a 
consequence of the human intervention in these areas. Carbayo 
et al. [5] also found lower abundance of land planarians in areas 
reforested with Pinus sp. than in native forests. Baptista and Leal- 
Zanchet [17] recorded lower abundance of flatworms in forests 
under some kind of anthropogenic influence. Furthermore, the 
composition of assemblages of planarians of the undisturbed sites 
of CAMB (Bl and B3) was similar to that of the conserved sites of 
RVSU (U1-U4). These six sites possibly provide greater 
availability of microhabitats to be colonized by land planarians 
[42]. In sites B2 and B4, the factors that would limit the presence 
of land planarians are the low complexity of the forest, the low 
availability of refuges in the soil, and the high incidence of sunlight 
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Figure 6. Non-metric multidimensional scaling (nmMDS) analysis based on the composition of land planarian assemblages of each 
sampling site. Sites with undisturbed condition (U1-U4, Bl and B3) were clustered at 54% similarity (Bray-Curtis coefficient), while disturbed sites 
B2 and B4 were segregated. 
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Table 5. Partitioning of gamma diversity in alpha and beta components for species richness (°D) and exponential Shannon 
diversity ('D) compared to null models (expected values) calculated between sites. 
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on the soil owing to the absence of forest strata and subsequent 
moisture loss and compaction of the soil [5,19,45]. 

The beta diversity between sites was higher than expected, 
indicating that the species turnover between sites contributed more 
to the total richness (gamma diversity) than the alpha diversity. 
The high beta diversity evidences that the sampling sites share few 
species [46] , probably as a consequence of the particular 
characteristics of each site of RVSU and CAMB, which harbor 
different species of land planarians and restrict many of them to 
certain areas. 

Our results suggest that land planarian assemblages are sensitive 
to the conservation degree of the rainforest. Our study also 
evidence that not only continuous areas with native forest (e.g., 
RVSU) but also fragmented landscapes in which the original 
rainforest is conserved as patches (e.g., Bl and B3 in CAMB) may 
harbor similar levels of land planarian diversity. Considering the 
high levels of beta diversity between sites and the low vagility of 
land flatworms, we suggest that conservation policies should 
promote the connectivity between fragmented areas. 
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